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Abstract 

The graph isomorphism problem is theoretically interesting and also 
has many practical applications. The best known classical algorithms 
for graph isomorphism all run in time super-polynomial in the size of the 
graph in the worst case. An interesting open problem is whether quantum 
computers can solve the graph isomorphism problem in polynomial time. 
In this paper, an algorithm is shown which can decide if two rooted trees 
are isomorphic in polynomial time. Although this problem is easy to 
solve efficiently on a classical computer, the techniques developed may 
be useful as a basis for quantum algorithms for deciding isomorphism 
of more interesting types of graphs. The related problem of quantum 
state symmetrization is also studied. A polynomial time algorithm for 
the problem of symmetrizing a set of orthonormal states over an arbitrary 
permutation group is shown. 

1 Introduction 

The problem of deciding if two graphs are isomorphic has many practical ap- 
plications such as searching for an unknown molecule in a chemical database 
[T] , verification of hierarchical circuits [2] and generating application specific in- 
struction sets [3] . It is not known to be solvable in polynomial time on a classical 
computer despite a great deal of effort to develop an efficient algorithm. On the 
other hand, its NP-completeness would imply the collapse of the polynomial 
hierarchy to the second level [I]. Moreover, the complexity class NP n coAM 
that contains graph isomorphism [5 is similar to NP n coNP which contains 
the decision version of factoring — a problem that is solvable in quantum poly- 
nomial time using Shor's algorithm [Bj. This suggests that graph isomorphism 
is of intermediate complexity and may be solvable in polynomial time on quan- 
tum computers. This paper takes a step towards a polynomial time quantum 
algorithm for graph isomorphism based on the state symmetrization approach 
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[7J by developing quantum algorithms for rooted tree isomorphism. By consid- 
ering all possible roots in one of the trees, it is also possible to efficiently decide 
if two unrooted trees are isomorphic. Although tree isomorphism can be de- 
cided in linear time on a classical computer [5] , the quantum algorithm for tree 
isomorphism shown in this paper relies on a different strategy which may be 
useful for solving graph isomorphism in cases that are not known to be solvable 
on classical computers in polynomial time. An algorithm is also presented for 
the related problem of symmetrizing a sequence of orthonormal states over an 
arbitrary permutation group. 

2 A quantum algorithm for tree isomorphism 

One approach to the graph isomorphism problem is to prepare a superposi- 
tion of all possible permutations of the adjacency matrix [TJ. For a graph 
G with automorphism group Aut(G) and adjacency matrix A, this state is 

\G) = y/^E^s. /Aut(G) where A n denotes the adjacency matrix ob- 
tained by applying it to A. Given two graphs G and H, the states \G) and \H) 
are equal if G and H are isomorphic and are orthogonal otherwise. It can then 
be determined if G and H are isomorphic using the swap test [5]. However, this 
is not sufficient to solve graph isomorphism since no efficient algorithm is known 
for preparing the state \G). As a first step towards developing a general quan- 
tum algorithm for graph isomorphism using this approach, a quantum algorithm 
for tree isomorphism based on state symmetrization will be presented. The tree 
isomorphism algorithm operates on qudits with four computational basis states 
rather than qubits. These states are |0), |1), \a) and The states |0) and 
|1) are used in most parts of the algorithm for representing binary numbers and 
permutations while the other computational basis states have special purposes. 
The state \a) is used to mark the beginning of a tree while |/3) marks the end 
of a tree as in the classical tree isomorphism algorithm. The algorithm works 
by preparing a state |T) for a given rooted tree T such that given two rooted 
trees T\ and Ta, if \T\) and |Ta) are of the same dimension then (TJT2) = I if 
and only if T x T 2 and (Tx\T 2 ) = if and only if T 1 ^T 2 . If |T X ) and |T 2 ) are 
not of the same dimension then it is required that T\ ^ T 2 . Such states can be 
prepared recursively as follows. Consider a rooted tree T. If T consists of a sin- 
gle node (for instance if T is a leaf node in a larger tree), then let |T) = \a, 0). 
Otherwise, let Ti, ... , T/. denote the subtrees rooted at the k children of the 
root node of T. For each 1 < i < k, let f/, denote the unitary matrix such that 
Ui |0) = |Tj). Note that the states |T) and \Tj) are not necessarily of the same 
dimension for Ti p= Tj as there are bijections between the number of qudits in 
|Tj), the dimension of |Tj) and the number of nodes in the tree Ti. Let bi denote 
the number of qudits in the state |Tj). Observe that (assuming the algorithm 
works properly on smaller trees) bi ^ bj or (Ti\Tj) = if Tj ^ Tj while bi = bj 
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and (Ti\Tj) = 1 if Tj = Tj. Now construct the state 

k 

<S> l T i>T, ® \°)au (!) 
t=i 

where each register Tj contains 6j qudits and each register Mi contains |~log 2 C\ 
qudits where I is the largest branching factor that occurs in the entire tree. Let 
Cj be the number of subtrees Tj such that rij = rij where rii is the number of 
nodes in Tj. Now execute algorithm [T] to obtain the state 

k 

(g) |T) Ti ® |m(i)) M< (2) 

i=l 

where m(i) is the number of subtrees Tj with j < i that are isomorphic to Tj. 
Now append a new set of registers initialized to |0) yielding the state 



Algorithm 1 Counting the number of subtrees Tj such that j < i and T = Tj 
for each 1 < i < k 



l: for i = 1, . . . , k do 

2: if Cj > 1 then 

3: for j = 1, . . . , i — 1 do 

4: if &, = 6j then 

5: Apply [/J to register Tj resulting in the state |0) if and only 
if Tj = Tj 

6: Add one to the contents of the register Mj conditional on the 
contents of register Tj being |0) 

7: Apply f/j to register Tj in order to restore the original state 

8: end if 

9: end for 

10: end if 

11: end for 



k k 

(g) |Tj) T . ® \m{i)) Mi (g) |0)^ (3) 
i=l i=l 

where each register A/j uses [log 2 qudits. XOR the contents of each register 
Mi into the register iVj obtaining the state 

k k 

(&\Ti) Ti ®\m(i)) Mi §§\m{i)) Ni (4) 

i=l i=l 

Letting |rj) = |T) (g) |m(z)) results in the state 

k k 

i=i »=i 
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where each register T[ contains registers Tj and Mj. Now, prepare the state \Sk) 
(the uniform superposition over all permutations in Sk) in another register to 
obtain 



fc fc 

S) \ t Jt> <S> l m «U ( 6 ) 
i=i j=i 



Apply a conditional permutation which permutes the order of the registers T[ 
resulting in the state 

fc fc 



J- y 



W ) T ,(g)|m(i))^ (7) 



Note that since the registers T[ now contain different numbers of qudits for differ- 
ent terms in the quantum superposition, special care must be taken here. Each 
7r G Sk is represented as a product <?fc • • • <?2 where each gi € {(1, i), . . . , (i — 1, i), t} 
and the permutation |vr) is stored as \gu) ® • • • ® \g%) where \gt) — \j) when 
9i = Cm)- The next step is to uncompute the contents of the S register. 
Each vector ® i=1 \T n -i^) T , corresponds to a unique permutation because each 

\n) = \Ti) ® \m(i)) is orthogonal to every \rj) — \Tj) <g> \m(j)) for i ^ j. Note 
that Ti might be isomorphic to Tj for some i ^ j but in this case m{i) ^ m(j) so 
that |rj) and \tj) will still be orthogonal. This was the reason for labeling each 
register Tj with the number of isomorphic subtrees in the sequence T\, . . . , Tj_i. 
Add another set of registers initialized to |0) obtaining the state 

■iE0 i°)p. Ms Kt-vo)^ (8) i m «U ( 8 ) 

^ ireSfe i=l i=l * t=l 

where each register Pj uses |~log 2 (fc + l)] qudits. Now XOR the contents of 
register Pi with \i) resulting in the state 

~m Y,®^P>®Ws®\^-H l ))T;®\ m ( t ))m (9) 

^ TreS k i=l i=l * i=l 

Apply a conditional permutation to reorder the registers P\,...,Pk according 
to the contents of register S obtaining the state 

TP E (8)l 7r_1 w)p^i 7r )s(8)l T -- i w)T;(8)i w W)^ ( 10 ) 

' ireSt i=l i=l i=l 

Now, note that each \n) is represented in subgroup tower form as a product 
5fc • • • ff2 where each gt e {(1, i), . . . , (i — l,i),t}. Then the contents of register S 
can be uncomputed using the procedure (to be described later) for converting the 
row representation ® i=1 |7T _1 (i)) p of n^ 1 to the subgroup tower representation 
\n) for 7r. This results in the state 
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ireSj, i=l i=l i=l 
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after discarding the contents of register S which is now |0). The next step is to 
uncompute each Pi register. To do this, first add a new register initialized to 
|0) yielding the state 

4y E l°)B0K" 1 (i)> F4 (g)k.- 1 ( i )> T/ (g)Ki)) % (12) 

' TrGSfc i=l i=l i=l 

where the register i? uses one qudit. The contents of the Pi registers are then 
uncomputed using algorithm [2] This results in the state 

after discarding the contents of the Pi and £? registers which are now all in the 
state |0). It is now necessary to uncompute the contents of the Ni registers. 



Algorithm 2 Uncomputing Pi for all 1 < i < k 
1: for i = 1, . . . , k do 
2: for j = 1, . . . , k do 

3: XOR the contents of register B with 1 conditional on register Tj 

containing bi qudits 
4: if a > 1 then 

5: Apply W to register Tj conditional on the contents of register B 

being |1) 

6: XOR the contents of register Pj with i conditional on the contents 

of register B being |1), the contents of register Tj being |0) and the contents 
of register Mj being equal to register Ni 

7: Apply Ui to register Tj conditional on the contents of register B 

being |1) to restore the register to its original state 

8: else 

9: XOR the contents of register Pj with i conditional on the contents 

of register B being |1) 
10: end if 

11: XOR the contents of register B with 1 conditional on register Tj 

containing bi qudits, restoring register B to the state |0) 
12: end for 
13: end for 



First add a new set of registers initialized to |0) yielding the state 

4fE<8) k-- 1 (o)i- i°)t» (8) m*)U ( 14 ) 

where each register T" uses bi qudits. Initialize the contents of each register 
T'i where a > 1 to \T) by applying Ui to each register T, where a > 1. This 
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results in the state 



K-Wt; ® I^V (g) M*)U (15) 

' irESk i=l <=1 £=1 

where |T/) = |T<) if c. ( > 1 and |T/) = |0) if c, = 1. The contents of the N. t 
registeres are then uncomputed by applying the inverse of algorithm [1] which 
was used to set each register Ni to \m(i)} in the first place. This yields the state 

A= E (8)I^-w)t/(8)i^)t» as) 

' TfGSk «=1 £=1 

after the Ni registers which are now all equal to |0) are discarded. Uncomputing 
the contents of the T" registers is accomplished by applying llj to each register 
T" where Cj > 1. Discarding the registers T" then results in the state 

k 

TP E ®V*-*®)t> ( 1? ) 

v ires k i=i 

Prepending a \a) and appending a to this state yields the final state 

|T) = \a) ® ( E ® I^-W)^ J ® ^ ( 18 ) 

' VireSfc £=1 / 

Given two rooted trees T-y and T2, the states |Tl) and |Tg) as defined in the above 
equation have the property that if \T\) and IT2) are of the same dimension then 
(Ti|T 2 ) = 1 if and only if Ti T 2 and (Ti|T 2 ) = if and only if T ^ T 2 as 
desired. 



3 Complexity analysis of the tree isomorphism 
algorithm 

In this section, it will be shown that the quantum tree isomorphism algorithm 
presented in section [2] is polynomial. To do this, the first step is to create a 
recurrence for the running time of the algorithm. This results in 

T(n) = sup {t ni ,...,n h } (19) 

EiU rii = n-l 
1 < fc < n- 1 
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where 



fc 

t ni ,...,n k = y ' T(jlj) 
i=l 

k i-1 

+ > i] I>* = n o\ ( 27 ><) + °( lo g^)) 

i=l j=l 

+ 0(fclog£) 

+ 0(k 2 \og£) 

+ Pi(S(n),k) 

+ 0(fclogfc) 

+ pi(k |Tog 2 (fe + 1)"| , fc) 

+ P2(k) 

k k 

+ E E ( 2 »»( n > fc ' & ') + & > !] ( 2T (^) + ( lo S^)) + [Ci = iMlog fc)) 
*=i i=i 
fe 

+ ]T[c 4 > l]T(ni) 

»=i 

fe i-1 

+ 1>* > = n A ( 27 ><) + °( lo g^)) 
»=i j=i 

fe 

+ > i]T(m) 
»=i 

+ 0(1) (20) 

and 

5(n) = 2n + riog 2 ^1 (n - 1) (21) 

is the number of qudits required for the state \T) where T is a tree containing n 
nodes. Each a denotes the multiplicity with which m occurs in the values rij as 
in the description of the algorithm. The functions pi describe the complexities 
of different subroutines of the algorithm and are bounded above by polynomials. 
The recurrence T(n) is bounded above by the recurrence 

T'(n)= ^ sup {/:, „. } (22) 

Ei=i n t =n-l 
\<k<n-l 

where 

k 

t' ni ,...,n k =J2(l + 2[ Cl >l}+ 4[cj > l]ci + 2fc[c, > 1]) T'(ni) + p(n, *) (23) 

i=l 
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and p is bounded above by some polynomial. It will now be shown that T'(n) 
is polynomial in n. Consider T"(n + 1). There exist k and iii with 1 < k < n and 

Eti m = n such that T'{n+l) = Eti (! + 2[ci > 1] + 4[ci > l]cj + 2fc[c 4 > 1]) T'(n 4 )+ 
p(n + 1, £). Let C be a subset of {1, . . . , k} such that for each rij, there exists 
j G C such that = rij and if i 7^ j g C then 7^ rij. Thus, for each distinct 
value of rij, C contains exactly one index j such that = rij. The analysis is 
performed using T'(n + 1) — T'{n). Consider the following three cases: 

• Suppose that nj = 1 for some 1 < j < k with j e C. One of the subtrees 
with only one node will now be removed to obtain a lower bound for T'(n). 



Let 




(24) 



Then a lower bound for T' 



(n) is 



T'(n) > 



£ (1 + 2[cJ > 1] + 4[c{ > l]cj + 2fc[cJ > l})T'( ni )+p(n,£) 




Therefore 
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T'(n + 1) - T'(n) < ^ (1 + 2[c; > 1] + 4[c, > l\a + 2k[ Cl > 1]) T'{m) +p(n + l,£) 

l<i<k 

- ]T {l + 2[c' i >l]+l[c\>l]c' i + 2k[c' i >l))T'{n i )-p{n,t) 
1 < i < k 

(26) 

= ^ c, (1 + 2[c 4 > 1] + 4[c, > l] Ci + 2k[ Cl > 1]) T'(ni) + p(n + lj) 

- 4 (1 + 2[c{ > 1] + 4[cJ > l]cj + 2*[cJ > 1]) r'(ni) - p(n, 

(27) 

< ^ c t (1 + 2[c 4 > 1] + 4[c, > l] Ci + 2k[c, > 1]) T'(m) +p(n + l,£) 

iec 

- Yl c i {l + 2[c i >l}+4[c i >l]c i + 2k[c i >l})T'(n i )-p(n,£) 
iec\{j} 

(28) 

=Cj (1 + 2[cj > 1] + 4[cj > l]cj + 2k[ Cj > 1]) T'irij) + p(n + 1,£) - p(n, £ 

(29) 

<n(3 + 6ra)T'(l) +p{n + !,(.)- p{n, i) (30) 

• Consider the case where Cj = 1 for some 1 < j < k. A lower bound for 
T'(n) is now obtained by removing one node from the subtree containing 
rij nodes. Let 

< = h~ l l rJ. (3D 

[n-i 

and let c- denote the multiplicity with which n- occurs in the values 
ni,...,n' k . Then a lower bound for T'(n) is 



T'(n) > ^(l + 2[c^ > 1] +4[4 > l]c' l + 2k[c' l > l])T'(n{)+p(n,^) (32) 
Thus, 



i=i 
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T'{n + 1) - T'{n) < ^ (1 + 2[c, > 1] + 4[c, > l]a + 2fc[c« > \])T\m) + p{n + l,i 



»=i 

k 



- ^ (1 + 2[4 > 1] + A[c' t > l]cj + 2fc[cJ > 1]) T'K) - p(n, £) 
»=i 

(33) 

it 

< ^ (1 + 2[ Ci > 1] + 4[ci > l]c< + 2fc[c, > 1]) T'(ni) + p(n + I, I) 



»=i 

fc 

-5^(l + 2[c< > l]+4[ci > l] Cl + 2fc[c, > l])T'(n;)-p(M) 

i=l 

(34) 

=T'(n,) - T'in'j) + p(n + 1,1) ~ p{n, t) (35) 

=T'{n 3 ) - T\n 3 - 1) + p(n + 1,£) - p(n, I) (36) 

< sup {T'(n')-T'(n' -l)}+p(n + l,e)-p(n,t) 
1 < ri < n 
n'efj 

(37) 

• For the final case, suppose that c, > 1 and rij > 1 for all 1 < i < k. 
Choose the smallest rij such that j e C. A lower bound for T'(n) will be 
obtained by reducing each m such that n, = % by one. Let 

n^h' 1 Hi=nj (38) 
[Hi rii^ rij 



Then a lower bound for T'(n) is 



T'(n) > (1 + 2[ci > 1] + 4[c, > l]ci + 2fc[c 4 > 1]) T'(nJ) +p(n, £) (39) 
»=i 

Hence, 
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T'{n + 1) - T'{n) < ^ (1 + 2[c, > 1] + 4[c, > l]a + 2fc[c 2 > 1]) T'(n;) + p{n + l,£) 



i=i 



- (! + 2[ci > 1] + 4[ci > l]c< + 2fc[c, > 1]) T'(n{) - p(n, 
»=i 

(40) 

= ^ Ci (l + 2[ Ci > l]+4[ Ci > l] Cl + 2fc[ Q > l])T'(ni)+p(n + l,0 



iec 



- 5^^(1 + 2^ > l]+4[ Ci > l] Cl + 2k[ Cl > l])T'«)-p(n,f) 

(41) 



= °i ( 3 + 4c * + 2fc ) T '(^) + P(n + 1, 



iec 



(3 + 4Q + 2fc) T'«) - p(n, 



(42) 



=cj (3 + 4c, + 2k) (T'irij) - T'(nJ)) + p(n + M) - p(n, £) 

(43) 

=cj (3 + 4c, + 2fc) (T'( nj ) - T'(n 3 - 1)) +p(n + 1,1) ~ p(n,l) 



(44) 



< 



sup {k('S + 6k)(T'( - ) -T'( 

k<k'<n-l 1 V VLfcJ/ V 

+ p(n + l,l)-p(n,l) 



'))} 



(45) 



A recurrence D(n) can now be denned in terms of maximum of the upper bounds 
obtained in each case so that 



Let 
where 



T'(n + 1) - T'(n) < D(n + 1) 
D(n + 1) = max{f(n),g(n),h(n)} + d(n + l,£) 

f(n) =n(3 + 6n)T'(l) 

<?(n) = sup {D(n')} 
1 < n 1 < n 
n' € N 

/i(n)= sup [fc(3 + 6fc)L»( ^ )} 

k<k'<n-l 1 VLfcJ/J 

fc' e Q 



(46) 
(47) 

(48) 
(49) 

(50) 
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and d(n + 1, £) is a polynomial upper bound for p(n + 1, t) — p{n, £). Also, note 
k > 2 is a value which depends on n. It will now be shown that D{n) < an v i q for 
all n and for all £ with appropriately chosen constants a, p and q. Let deg n d and 
deg^ d denote the degrees of d in n and £ respectively. Let p = max{dcg n d+1, 6} 
and q = deg e d. Since (n + l) p = J]fe=o (S nfe ' ( n + X ) P ~ nP > pn"' 1 ■ This 
implies that by choosing a constant a > which depends only on the coefficients 
of d(n,£) and the constants D(l) and T'(l), the following properties can be 
satisfied: 

• D(l) < a 

• f( n ) < an p ^ 9 for all n and for all £ 

• d(n + !,£)< a{n + l) p £ q - an p i q for all n and for all I 

It can then be shown that D(n) < an p £ q for all n and for all £ by induction. 
Note that the basis case is satisfied by the choice of a. For the inductive case, 
assume that D(n') < an' p £ q for n' < n. Now consider D(n + 1). By choice of 
a, f{n) < anP£ q so 

f(n)+d(n+l,£) <a(n+lf£ q (51) 



By assumption, D(n') < an' p £ q so for 1 < n' < n, D(n') < an p £ q . Thus, 
g(n) + d(n + l,£)<a(n+l) p £ q 



(52) 



Recall that in the expression for /i(n), k > 2. Therefore, for k < k' < n— 1 with 

k' G Q, 



fc(3 + 6fc)L> ( 



<afc(3 + 6A;) - 
Lfc'_ 

<afc(3 + 6fc) (-J 
<an p £ q 



Thus, /i(n) < an^« so 

Consequently, 
which proves that 



h(n) + d(n + l,£)<a(n+l) p £ q 

D(n+l) < a(n+l) p £ q 
D(n) < an p £ q 



(53) 

(54) 

(55) 
(56) 

(57) 

(58) 
(59) 
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for all n and for all I. Hence, 



T(n) <T\n) (60) 
n 

= ^2(T'{i)-T'(i-l)) (61) 

i=l 

■n 

<Y,m (62) 

i=l 

n 

<a^« ^ z p (63) 

n 

<a^V^ (64) 

=an p+1 ^ (65) 
<an p+q+1 (66) 
=0(n p+9+1 ) (67) 

so T"(n) is polynomial in n. 

4 Subroutines for the tree isomorphism algo- 
rithm 

To finish describing the algorithm, all that remains is to show how to implement 
the three subroutines it utilizes. These are preparing the uniform superposition 
over all permutations, applying a conditional permutation and converting the 
row form 0*^ \TT-\i)) of the inverse of a permutation 7r £ Sk to the subgroup 
tower form \n) if n. 

4.1 Preparing the uniform superposition over all permu- 
tations 

This state can be prepared for Sk using 0(k 2 log k) basic operations using gener- 
ator state based algorithms for initialization [101 HU [H] • It is actually possibly 
to do this more efficiently but this is not necessary for the purposes of this 
algorithm as it will result in the same complexity. 

4.2 Implementing conditional permutations 

Since this operation is a bijection which takes an input of fixed length that 
can be implemented in polynomial time on a classical computer it can also be 
implemented in polynomial time on a quantum computer. 
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4.3 Converting permutations from row form to subgroup 
tower form 

This section will show how to convert a permutation n € Sk in row form 

k 

®k(i)> fl (68) 

z=l 

to the subgroup tower form 1 7r) . First, observe that (7r(fc),fc)7r 6 Sfc-i so the 
equation 

7T = (7r(fc),fc) ((7r(fc),fc)7r) (69) 

allows 7r to be decomposed into the product of one of the transpositions in 
the set {(1, k), . . . ,(k — 1, k), l] and an element of Sk-i- The conversion proce- 
dure works by repeating this decomposition. Append a second set of registers 
initialized to |0) resulting in the state 

k 2 

®k(i)>« ®|0> G< (70) 

i—l i—k 

where the register Gi contains |~log 2 (« + 1)] qudits. Now execute algorithm [3J 
This results in the state 

k 2 

0N)p 4 (8)|5fe) Gi (71) 

i— 1 i—k 

where 7r = ■ ■ ■ gi is the subgroup tower form of n. Discarding the registers 

Algorithm 3 Converting a permutation in row form to subgroup tower form 
1: for i — k, ... ,2 do 

2: Copy the contents of register Pi to register Gi 
3: for j = 1, . . . , k do 

4: Swap the contents of registers Pi and Pj conditional on j being equal 

to the contents of register Gi 
5: end for 
6: end for 



Pi yields 

2 

K) = (g) l<7fc) G< (72) 

i—k 

which is the subgroup tower form of n. Now the algorithm actually requires a 
procedure that will convert the row form 

k 

(E)K _1 w)p. ( 73 ) 



14 



of 7r 1 to the subgroup tower form | vr) of ir. To do this it suffices to show how 
to invert the row form 

k 

®k«>P 4 (74) 

of 7T. This can be done by first appending new registers each in the state |0). 
This results in the state 

k k 

®k(i)>* ®|0) o , (75) 

i=l i=\ 

where each register Qi contains |~log 2 (fc + 1)] qudits. Then apply algorithm 2J 
This results in the state 

k k 

<g)\K(i)) Pi <g)\K- X (i)) Qi (76) 
i=l i=l 

Applying algorithm 2] again with the roles of the registers Pi and Qi swapped 
Algorithm 4 Computing the inverse of the row representation of tt 



for i = 1, . . . , k do 



for j = 1 , . . . , k do 

XOR the contents of register Qj with \i) conditional on register Pi 
being equal to \j) 
4: end for 
5: end for 



uncomputes the contents of the Pi registers. This yields the state 



k k 

^)pM\^\i))o i ( 77 ) 



i=i i=i 
which results in the row form 

k 



l^ 1 (0) o . (78) 



after discarding the registers Pj. The cost of this procedure is 0(k 2 logfc) basic 
operations. 

5 Quantum algorithms for state symmetrization 

In this section, the problem of quantum state symmetrization which is closely 
related to graph isomorphism will be studied. Let , . . . , \ip n ) be a sequence 
of ortho normal states and let G be a subgroup of S n . The problem is to 
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prepare the state — j= ^2 nGG (Vv- 1 ^))- Consider a graph G with ad- 

jacency matrix A. Suppose that it was possible to efficiently prepare the state 

\G) = v^i^^ 

/Aut(G) \^*) where Aut(G) denotes the automorphism 
group of G and A* denotes the adjacency matrix obtained by applying ir to 
A. Then given two graphs G and H, one could create the states |G) and \H) 
and these states would be the same if G and H were isomorphic and orthogonal 
otherwise. The swap test [9] could then be used to determine if G and H were 
isomorphic [7j. The only difference between the state symmetrization problem 
and the state preparation approach to graph isomorphism is that in the former, 
symmetrization is performed over a sequence of orthogonal states whereas in 
the latter symmetrization is performed on the adjacency matrix. Moreover, the 
quantum algorithm for tree isomorphism shown in section [2] works by repeatedly 
symmetrizing over Sk- An algorithm will now be presented for performing state 
symmetrization over an arbitrary permutation group G on n elements. 



5.1 A quantum algorithm for state symmetrization over a 
permutation group G on n elements 

Let G be an arbitrary permutation group on n elements for which a gener- 
ating set K is known. It will be shown how to efficiently prepare the state 
-7= J2tt£G ®i=i | V'tt— 1 (i) ) where \tpi) , . . . , is a sequence of orthonormal 

states as before. The algorithm works using a classical group theoretic algo- 
rithm as a subroutine. Consider the tower of subgroups 

G (0) > G (l) > > Q{n) ( 7Q ) 

where G^ = Gi j with G\ i denoting the pointwise stabilizer of the set 
{1, . . . , i} in G. Note that G^ is the trivial group. Now let Ui be a left transver- 
sal of GW in G^ l ~ l > — that is, a set containing exactly one representative of 
each left coset of G^ in G^ -1 '. Then observe that each of the permutations in 
G may be uniquely expressed as a product of the form 

9i ■ ■ ■ 9n (80) 

where each g. L £ Ui . This can be thought of a generalized version of the subgroup 
tower form utilized in the quantum algorithm for tree isomorphism. Now, there 
is a classical algorithm [13j [14] which can find the left transversals Ui given 
the generating set K for G using 0(|if| n 2 + n 6 ) operations. Moreover, it can 
be shown that the index of G^ in G 1 - 1-1 - 1 is at most n — i + 1. This allows 
symmetrization over G to be performed efficiently as follows. First, create the 
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state 



1 n 1 

M i=l V 1^*1 ftSC/; 

1 " 

2 ®lft) Gl (82) 



g% e £/* 

where 1 7r) denotes the permutation 7r represented as a product of the form shown 
in equation (|80[) . Note that this can be done efficiently since it is easy to prepare 
each of the states — ^= j E 9 e(7 l^*) an< ^ * ne desired state is simply their tensor 
product. The representation used for each permutation 7r will now be changed 
to row form. Add a new set of registers yielding the state 

n n 

^ ®ISi>G 4 ®l*>J>, (83) 



l G l 

n = 9f-9n 
9i G Ui 

where each register Pi is initialized to \i) and uses |~log 2 (n + 1)] qubits. The next 
step is to compute the row form of each permutation g% ■ ■ ■ g n in the registers 
Pi . This is done using algorithm [5] and results in the state 

- n n 

E ®\9i) at &\*®)p i (84) 



\G\ 

7T = 91 ■ ■ ■ 9n 1 

9i € Ui 

It is now necessary to uncompute the contents of the Gi registers. Now, given 
Algorithm 5 Computing the row form of each permutation 



for i — 1, . . . , n do 



for gi G Ui do 

Multiply the permutation represented in row form by the registers 
Pi , . . . , P n on the right by gi conditional on the contents of register Gi being 
9i 

4: end for 
5: end for 



the left transversals Ui, there is an efficient classical algorithm [14] which takes 
a permutation ir e G and computes values gi E Ui such that tt = g\ ■ ■ ■ g n . 
The left transversals Ui can be regarded as fixed so that the only input of the 
algorithm is the permutation tt. The algorithm therefore takes an input of 
fixed size so it can be converted to a classical circuit of polynomial size which 
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computes the values <?;. This circuit can therefore be utilized to uncompute the 
contents of each Gi register which yields the state 

1 " 

E (85) 



7T = 91 ■ ■ ■ 9n 

gt € Ui 

after the Gi registers are discarded. This can be rewritten as 



n 

1 E^W 1 »p, (86) 



Now, append the orthonormal states \ipi) , • ■ • , \ipn)- This results in the state 

1 n n 

pE<8>W'))n0WB i (87) 



ttSG i=l i=l 

Applying a permutation to the order of the P; registers conditional on the 
permutation represented by the registers Pi results in 

n n 

= £(g)M0}p s (g)K- lw > Bl (88) 

\ u \ TrgG i=l i=l 

Next, use algorithmic to invert each permutation tt represented by the registers 
Pi. This yields the state 

^ n n 

^fE®!^))^®!^-^ (89) 

Algorithm |5] is then applied to uncompute the contents of the P, registers. The 
result is 

Algorithm 6 Uncomputing Pj for all 1 < i < n 
l: for i = 1, . . . , n do 
2: for j = 1, . . . , n do 

3: XOR the contents of register Pj with i conditional on the contents 

of register Bj being \tpi) 
4: end for 
5: end for 



^ n n 



ttSG i=l i=l 

which yields the desired state 



1 ™ 

I I 7T£G 4=1 



(90) 



(91) 



after the P,; registers are discarded. 



18 



6 Conclusion 



This paper presented a new quantum algorithm which is capable of deciding 
tree isomorphism in polynomial time; although there is no difficulty in decid- 
ing tree isomorphism in polynomial time on classical computers, this quantum 
algorithm relies on new techniques which may be useful for more general quan- 
tum algorithms. The state symmetrization problem was also discussed and a 
quantum algorithm was shown which can be used to symmetrize a sequence of 
n orthonormal states over any subgroup of S n . 
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